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«TITLE MTHSHTANH : W Floating Hyperbolic Tangent routine 
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ALL RIGHTS RESERVED 
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N ACCO aoe CE WITH THE TERMS OF SUCH LICENSE AND WITH THE 
INCLUSION OF THE ABOVE yp lsh NOTICE. THIS SOFTWARE OR ANY OTHER 


5 
® 
® 
® 
® 
® 
® 
® 
* 
R Y PROVI » 
OTHER PERSON. NO TITLE TO AND OWNERSHIP OF THE SOFTWARE IS HEREBY * 
am 
® 
® 
® 
® 
® 
& 
® 
® 
® 
® 


* 
.* 
ft 
* 
* 
* 
: 
:f 
:& 
:* TRANSFERRED. 
. 
:f 
+f 
+f 
st 
ft 
te 
.f 
% 


o 
za 
a 
«< 
_— 
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fACILITY: MATH LIBRARY 
"ABSTRACT: 


MTHSHTANH is a function which returns the H floating hyperbolic tangent 
of its H floating argument. The call is standard call-by-reference. 


VERSION: 1 
HISTORY: 
AUTHOR: 
John A. Wheeler, 20-Oct-1979: Version 1 


MODIFIED BY: 
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-SBTTL HISTORY ; Detailed Current Edit History 


MTHSHTANH ; HF 
1-006 HISTO 


Edit History for Version 1 of MTHSHTANH 


1-00) - Adapted from MTHS$GTANH version 1-002. JAW 20-Oct-1979 
1-002 = Change constant 16.0 to 4 to correct inaccuracy. The 
value of X above which 1.0 is the best machine a proximation 
to HTANH(X) is about 39.86. The next higher num gr that can 
be represented as a short Literal is 40. JAW 19-Sep-80 
Eliminated syave; i short literals. RNH 15-Oct-81 
Changed #40 to S*#40 to ¢tecumyent gssenb\ er problem RNH 23-0ct-81 
Use general mode addressing. SBL_ 30-Nov-1981 
Changed the constant 2*-59 8 Au5 te correct inaccuracy. For 
values of iX! between 2°-59 and 2°=-57 no noticeable loss of 
sront scence was noticed by ghe assumption that HTANH(x)=x starting 
at iXi<=2*=59 instead of 2°-57, but a loss of performance was felt 
do to unnecessary eae of HSINH(X)/HCOSH(X). ALL soproge late 
references to 2*=59 have been changed to 2°-57. JCW 10-Jan-1983 
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«SBTTL DECLARATIONS 


“oO 


INCLUDE FILES: 


EXTERNAL SYMBOLS: MTHSJACKET_HDLR 


SABL GBL 

KTRN  MTHSHCOSH 
SEXTRN MTHSHEXP_R6 

; EQUATED SYMBOLS: 


; PSECT DECLARATIONS: 


. LONG 
H_2_POWER_M57: 
. CONG 


AUF WIN OC OONAUS WN 9 ODNAU EWN OO OONOULS WO 


SOGDOOOO0O000 000000000009 09000000090909 INI IIIS 


5 a SS a 
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MTHRTL.SRCJMTHHTANH.MAR; 1 


; Declarative Part of Module 


-EXTRN on all symbols 


; Offset for first argument 


PIC,SHR,LONG, EXE ,NOWRT 


Program section for math routines 


; 0.25 


208-57 
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ating Point TANH 1:25: MTHRTL.S RCIM THHTANH. MAR; 1 (4) 
~SBTTL MTHSHTANH = H Floating Point TANH 


PEED D GMs ££ EE BEE EWI WIP POPUNIPININPINININ 2 2 PS OO Os "OO “2 =x 
FWD CO OONAU EWR 9 OO NAW WN 9 OO NAUE WIN (OOO NOU EWN — CO OONOAOUES WN OwWO 


artist 


dy 

; FUNCTIONAL DESCRIPTION: 

: HTANH - H floating point hyperbolic tangent function 
; HTANH(X) is computed as: 


If Xi <= 20-57 e the. HTANH(X) = 

If 2**-57 < IX} <= then HTANHGX) = HSINH(X)/HCOSH(X). 

If 0.25 < !x! < 40.0, oe HTANH(X) = (HEXP(2*X) = 1) / (HEXP(2*x) + 1) 
If 40.0 <= iXi, then HTANH(X) = sign(x) * 1 


2 CALLING SEQUENCE: 

htanh.wh.v = MTHSHTANH(x.rh.r) 

-or- 

CALL MTHSHTANH(htanh.wh.r, x.rh.r) 

Because an H-floating result cannot be expressed in 64 bits, it is 


returned as the first argument, with the input parameter displaced 
to the second argument, in accordance with the Procedure Calling 


Standard. 

: INPUT PARAMETERS: 

$4 44494 Seog : pet ine Longuor¢ multiplier 

s ; Contents of x is the argument 

: IMPLICIT INPUTS: none 
; OUTPUT PARAMETERS: 
; VALUE: H floating hyperbolic tangent of the argument 

00000004 : htanh = 1 * LONG ; htanh is the result 
IMPLICIT OUTPUTS: none 
; COMPLETION CODES: none 
: SIDE EFFECTS: none 
: NOTE: This procedure disables floating point underflow, enables integer 
: overflow. 

4OFC ENTRY MTHSHTANH, “M<IV, R2, R3, R4, R5, RE, 


R7> 
: Standard cat l= by-reference entry 
Disable DV (and FU), enable IV 
MTHSFLAG_JACKET : Flag that this is a jacket procedure in 
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a kk kk dk dh dd 
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: Qee-57 < 1X! < 40.0 


S eeeS7 < IX! <= 0.25 


MOVAB G*MTHSSJACKET_HND, (FP) 


; 2ee-57 < IX! 


CMPH RO 


BGEQ GEG_T0_40.0 


CMPH 


RO, H_0.25 
BLEQ LEQ_10_0.25 


0.25 < ix! < 40.0 


ADDH3 = ax 
JSB G* 
ADDH3S = #1 
SUBH2 = #1 
DIVH2 4 


RG; R 
BRB RETURN 


LEQ_T0_0.25: 
CALLG (AP) 


V me 
CALLG (AP), G6*MTHSHSINH 
MOVO ay (AP), RO 
R4, RO 


BRB RETURN 


GEQ_10_40.0: 
MOVH 


#1, RO 
TSTH ax CAP) 


41:24 VAX/VMS Ma 
$h:8% SR 


x 
CMPH = RO 
OUT. 
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cro 
MTHRTL.SRCJMTHHTANH.MAR; 1 (4) 


handler 


case of an error in routine 


set handler address to jacket | 
If an error, convert signal to user PC | 


and resigna 
O/R3 = X = @x(AP) 
/R3 = ix! 
Compare iX: with tae 
Branch if iXi <= 2*#-57 
; Compare iX: with 40.0 
: Branch if iX: >= 40.0 


Compare :X: with 0.25 
Branch if ix: <= 0.25 


RO/R3 = 2*X 

RO/R3 = HEXP(2*x) 

R4/R7 = HEXP(2*X) + 1 

RO/R3 = HEXP(2*X) = 1 

RO/R3 = (HEXP(2*X) = 1) / (HEXP(2*xX) + 1) 
Store result and return 

@y(AP) = HCOSH(X) 

R4/R7 = HCOSH(X) 

@y(AP) = HSINH(X) 

RO/R3 = HSINH(X) | 
RO/R3 = HSINH(X) / HCOSH(X) 

Store result and return | 
RO/R3 = 1.0 


Test the sign of X 


; 
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1-006 ATHSHTANN 2 WPloocing Polee TANR or eee tba f1:38i8¢ EMTART eaedarictana.mars1 22° (fy | 
4 18 y 17 pore RETURN 3 eronen if X >= 0 | 
50 8 72FD 8 18 EGH RO, R ; RO/R aa 
04 BC 7DFD F 219 RETURN: MOVO RO, @htanh(AP) ; Store result in first argument | 
04 3¢ 9 RET 3; Return to caller 
9 ; 
9 ¢ ; IMi <= 2e0=57 
38 55e | 
04 BC 08 BC 7DFD 9 § OUT_X: povo @x(AP), ahtanh(AP) ; Store result in first argument 
04 ge ; Return to caller 
9C § «END 


ee 
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Symbol table 7:34 (CMTHRTL.S RIM THHTANH.MAR; 1 (4) 


GEQ_10_40.0 0 p08 1R 01 
HTANH = 00 4 
“ 000 8 RB 1 
4 TOO. M57 4 10 R 1 
LEQ-10_0.25 8 88 $8 ge 1 
G = 00000004 
MTHSSJACKET_HND seeeereer Xx (1 
THSHCOSH eererere § xX 8 
MTHSHEXP_R6 rererere x 0 
eererere x 00 
MTHSHTANH 00000020 RG 01 
444 5 01 
QO008F 01 
= 00000008 
= 00000004 
pow maarece em eee ee fe 
: 3 _Psect synopsis ; 
PSECT name Allocation PSECT No. Attributes 
» OS ; 00000000 ( 0.) 00 ¢ O.) NOPIC USR CON ABS Len _ NOEXE NORD NOWRT NOVEC BYTE 
“MTHSCODE 0000009¢ (¢ 156.) 01 ¢ 1.) PIC USR CON REL LCL SHR EXE RD NOWRT NOVEC LONG 
$ccseecucsssunenesunassen} 
! ; Performance indicators ! 
Page faults CPU Time Elapsed Time 
Initialization 30 00:00:00.07 00:00:01.03 
roanens processing 117 00:00:00.64 00:00:04.22 
88 00:00:00.70 00:00:05.40 
Syabol table sort 0 BDO On Oe 00:00:00.00 
ass 57 sb S44 00:00:02.68 
oa Ee table output 2 00:00:00.03 00:00:00.04 
Psect synopsis output 2 00:00:00.02 00:00:00.06 
Cross-reference output 0 00:00:00.00 bs Bae Be 
Assembler run totals 298 00:00:02.10 00:00:13.4 


The working get Limit yes 750 pages. 
3485 bytes (7 pages) of virtual memory were used to buffer the intermediate code. 
oe were 10 pages of symbol table space allocated to hold 15 non-local and 0 local symbols. 
source Lines were read in Pass 1, producing 11 object records in Pass 2. 
1 page of virtual memory was used to define 1 macro. 


pewoeeceeseecescocceseseoces 


Macro Library name Macros defined 


-$255$DUA28: CSYSLIBISTARLET.MLB; 2 0 
0 GETS were required to define 0 macros. 
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VAX=11 Macro Run Statistics 
There were no errors, warnings or information messages. 


SEP=1984 


MACRO/ENABLE=SUPPRESSION/DI SABLE=(GLOBAL , TRACEBACK) /LIS=LIS$:MTHHTANH/ OBJ=OBJ$:MTHHTANH MSRC$:MTHJACKET/UPDATE=(ENHS$: MTHJACKET) #MSRC 
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